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Abstract 

Triple-quark elastic scattering amplitudes from perturbative QCD are first cal- 
culated and then used in a transport equation to study the thermalization of quark 
matter. By examining momentum isotropy to which the transport equation leads, 
we can determine thermalization time and offer an initial thermal quark distribu- 
tion function. With an anisotropic initial quark distribution, which is relevant to 
quark matter initially created in a central Au-Au collision at ^snn = 200 GeV, 
the transport equation gives a time of the order of 1.8 fm/c for quark matter itself 
to thermalize by the triple-quark scatterings. 

PACS codes: 24.85. +p; 12.38.Mh; 25.75.Dw; 25.75.Gz 
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1. Introduction 



Hydrodynamic calculations [1-7] that explain the elliptic flow coefficient V2 [8,9] 
at p± < 2 GeV/c, lead to the conclusion that quark-gluon matter created in Au-Au 
collisions at the Relativistic Heavy Ion Collider (RHIC) thermalizes rapidly. The usual 
parton scatterings give a long thermalization time of quark-gluon matter [10]. However, 
triple-gluon elastic scatterings in gluon matter lead to a short thermalization time of the 
order of 0.45 fm/c [11]. Such triple-gluon elastic scatterings can occur in view of the 
high gluon number density that has been achieved in Au-Au collisions [12-14]. However, 
since quark-gluon matter includes quarks and measured quantities rely on the presence 
of quarks, it is important to investigate the thermalization of quark matter based on 
initial anisotropic quark distributions. At present, the thermalization mechanism of quark 
matter is unknown. In this letter, we study the thermalization based on two- and three- 
quark elastic scattering processes. Until now, the squared amplitudes for the triple-quark 
scatterings have not been studied and must be calculated in perturbative QCD. Therefore, 
this work comprises two new tasks: one is to derive the squared amplitudes and the other 
is to study the thermalization. In order to keep the physics relating to the triple-quark 
elastic scatterings transparent, and due to the complications of calculating the squared 
amplitudes and numerically solving a transport equation, we exclude here the quark 
scatterings with gluons and/or antiquarks. 



2. Quark-quark and triple-quark scatterings 



The evolution of quark matter relies on the scattering of two quarks with identical 
or different flavors. Expressed in terms of the Mandelstam variables (s,t,u), the spin- 
and color-averaged squared amplitudes of order a 2 in perturbative QCD [15,16] are 
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for the scattering of up and down quarks. Here, g s is the quark-gluon coupling constant 
and gg = A-Kd s . 

The diagrams for elastic scattering processes of three quarks at order a* are drawn 
in Figs. 1 and 2. The three-quark scattering processes in Fig. 1 and Fig. 2 are separated 
into two distinct classes of diagrams: 

(1) A_, A_ 45 , A_ 56 , A_ 64 , A_4( 56 ), A_(45) 6 of Fig. 1 

involving two successive gluon exchanges and, 

(2) A*, A*45, A*56, A*64, A*4(5 6 ), A*(45) 6 in Fig. 2, 
involving a triple gluon coupling. 

When two final quarks have the same flavor, an exchange of the two quarks generates 
a different diagram from those diagrams A_ in Fig. 1 or A* in Fig. 2. When three quarks 
have the same flavor, exchanges of the three quarks lead to the diagrams A_4(56) and 
A_(45)6 of Fig. 1 and A* 4(56 ) and A*( 45)6 of Fig. 2. 
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A quark four-momentum is labeled as p\ = (E^pi) in the process q(pi) + q(p 2 ) + 
q(P3) ~~ * Q.(Pa) + ^(Ps) + q(P6)- The six momentum variables are assigned to quarks in six 
orderings. For instance, even though the momenta pi, P2 and p 3 may be assigned to the 
initial quarks from left to right of a diagram, while p 4 , p$ and p e are assigned to the final 
quarks from left to right, the initial quarks may also take the order of (pi, p 3 , p 2 ) and 
correspondingly the final quarks have momenta in the order of (p 4 , p e , p 5 ). If these two 
orderings are denoted as 

( Pa P5 Pe \ ( Pa Pe P5 \ 
V Pi P2 Pz ) ' V Pi Ps Pi)' 

the other four orderings are written as 

( P5 Pa Pe \ ( P5 Pe Pa\ ( Pe Pa P5 \ ( P6 Ps Pa \ 

\P2 Pi P3 J ' V P2 P3 Pi J ' V P3 Pi P2 J ' V ^3 P2 Pi)' 

The six orderings applied to the six diagrams in Fig. 1 generate 36 different diagrams with 
the corresponding momentum labels. However, the six diagrams in Fig. 2 are unchanged 
upon using the six orderings. Therefore, 42 different diagrams with various momentum 
labels are needed to take into account the squared amplitudes. 
Nine Lorentz-invariant variables are defined as 

Sl2 = (Pi + P2) 2 , S 23 = (pa + p 3 ) 2 , S31 = (p 3 + Pi) 2 

U12 = U15 = (Pi - PbY, W13 = M ie = (Pi - Pe) 2 

U21 = U2A = (P2 - Pa) 2 , u 2 3 = u 2e = (P2 - Pe) 2 

u 3 i = u u = (p 3 - p 4 ) 2 , U 32 = M35 = (p 3 - p 5 ) 2 

The squared amplitudes and interference terms are expressed in terms of these indepen- 
dent variables for the three-quark to three-quark scatterings. Our Fortran code includes 
the squared amplitudes summed over spin and color states of the three final quarks and 
averaged over the spin and color states of the three initial quarks. All interference terms 
of amplitudes of different diagrams are also calculated. The two interference terms of a 
diagram in Fig. 1 and a diagram in Fig. 2 cancel each other because the terms are pure 
imaginary. For example, 

Mt_ 45 M A *A(5e) + -Ma-45-M 1,4(56) = °- ( 3 ) 

We present the individually squared amplitudes and some interference terms of the 
second class of diagrams in Fig. 2, which are not zero, as 

I M A , | 2 = 3A/(w 2 lW 2 2 w 2 3 ), 
3 

3 _ 

\M AtA5 | 2 = SD/(uj 2 u 2 21 w 2 3 ), 
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Ma^MX^ = --F/(u 12 u 21 u 13 u 31 w 2 w 3 ), 

I M A ^ | 2 = SG/iulzulzwl), 
3 

M A* x MX t64 = --H / (u 23 u 32 u l3 u 3 iw l w 2 ), 
I M A , m | 2 = 3I/(u 2 13 u 2 31 w 2 2 ), 

I -^A^ (8a) | 2 = 3J/(Wi3«21«32)> 

3 

I ■ M A, (45)6 | 2 = 3L/(u? 2 u|,t4), 

where = Si 2 + S3i+«12+Wi3, ^2 = Si2 + s 23 +M 2 i+ M 23 and w 3 = S3i + S23+""3i+M32- A is 
a function of the nine Lorentz-invariant variables, A — A(ni, n 2 , n 3 , n 4 , n 5 , n 6 , n 7 , n 8 , n 9 ) 
Xs i2 s 23 s 3i M i2 M il M 2i M 23 M 3i M 32 where the sum is over rii, n 2 , n 3 , n 4 , n 5 , n 6 , n 7 , n 8 and 
n 9 , which run from to 4 and are restricted to give n\ + n 2 + n 3 + n 4 + n 5 + n 6 + n 7 + 
n 8 + n 9 = 4. The non-zero A's are tabulated to form the column under A in Table 1. 
B,C,D,E,F,G,H,I,J,K and Z are similarly defined and B,C,D,E,F,G,H,I,J,K 
and L can be found in Table 1. The above interference terms give the other interference 
terms, for example, MX,Ma, 4{56) = MaMa^y 

3. Thermalization of quark matter 

To study the thermalization of quark matter, these three-quark to three-quark scat- 
terings are incorporated into a Boltzmann-type transport equation 
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x [/i/ 2 / s (l - / 4 )(1 - / 5 )(1 - /e) - /4/ 5 /e(l - /i)(l - / 2 )(1 " / 3 )] 

(4) 

where the degeneracy factor gQ = 6 and the velocity of a massless quark vi = 1. The 
distribution function /; depends on the position rl, the momentum pi and the time t. The 
distribution function f\ in this equation gives the up-quark distribution in quark matter. 
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We assume that quark matter studied here has the same down-quark distribution as the 
up-quark distribution. The first term on the right-hand side of the above equation is 
for the 2-quark to 2-quark scatterings. The second term is a new term representing the 
3-quark to 3-quark scatterings and involving six quark distributions. 

The squared amplitude | -M U uu^uuu | 2 is for the scatterings of three identical quarks. 
It equals the sum of the individually squared amplitudes for the 42 diagrams discussed in 
the last section and the amplitudes corresponding to the interference terms between these 
diagrams. Both of the squared amplitudes, | -M uu d^ U ud | 2 an d | -A/f u d u ^udu \ 2 , are for the 
scatterings of two up-quarks and a down-quark. While the former receives contributions 
only from the diagrams A_, A_ 45 , A* and A* 45 , the latter has contributions only from 
the diagrams A_, A_ 64 , A* and A* 64 . The squared amplitude | -M u dd^udd | 2 is for the 
scatterings of one up-quark and two down-quarks. Only the diagrams A_, A_ 56 , A* and 
A* 56 contribute to the squared amplitude. 



4. Numerical results 



Anisotropic parton momentum distributions are formed in initial Au-Au collisions. 
Such anisotropy can be eliminated by elastic scatterings among partons. Starting from 
the time t ini , when anisotropic quark matter is formed and ending at the time t iso , when 
local momentum isotropy is established, the transport equation is solved. The squared 
amplitudes are calculated at a s = 0.6. In the calculations, a divergence of the squared 
amplitudes is encountered. The divergence is removed while gluon propagators are regu- 
larized by a screening mass which is evaluated from the distribution function by a formula 
used in Refs. [17-19]. We arbitrarily add 100 MeV to or subtract 100 MeV from the screen- 
ing mass, a change of thermalization time by about 40% is obtained. Any two runs of 
Fortran codes for solving the transport equation do not give rise to a difference of screen- 
ing masses more than 100 MeV. Therefore, the uncertainty in the thermalization time 
due to the screening mass is less than 40%. 

Since the initial gluon rapidity density accounting for the RHIC data [12,13] is about 
1000, larger than 200 given by HIJING [22], the initial quark rapidity density must also 
be similarly underestimated by HIJING simulation. For a central Au-Au collision at 
v /sjvtv — 200 GeV, an initial quark momentum distribution can be, for example, four 
times [23] as large as that obtained from HIJING simulation [24], 

^ = ff2 J. I | f^uZ n ^ e-|Pl/(0- 9 cos h ( y ))-(| P 1tan h ( y ))V8g- (y2 _ 2) (5) 

7vR A Y(\ p I / cosh(y) + 0.3) 

where | p | is in GeV, the nuclear radius Ra = 6.4 fm, t ini = 0.2 fm/c and the rapidity y 
is less than Y which approximately equals 5. 6(x) is the step function which is zero for 
x < or 1 for x > 0. Anisotropy of the distribution is shown in Fig. 3 by the differences of 
dotted, dashed and dot-dashed curves, which represent the quark momentum distributions 
at the three angles 9 = 0°,45°,90° relative to one incoming gold beam direction. The 
quark momentum distribution in the longitudinal direction is considerably larger than 
one in the transverse direction when momentum departs from 0. 

Numerical simulations of the 2-^2 and 3^3 scatterings in Fortran codes depend 
on quark positions. Starting from the cross section for the 2-quark to 2-quark scattering 
o" 2 ^2 which is calculated from the squared amplitudes in Eqs. (l)-(2), a scattering of two 
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quarks is considered to occur when the closest distance of the two quarks is less than an 
interaction range of yavV^". Regarding to a three-quark scattering, the three quarks 
must be in a sphere which center is at the center-of-mass of the three quarks and which 
radius r\ ls is determined by [11] 

nThs 3\J (2vr) 3 2£ 4 (2tt) 3 2£ 5 (2tt) 3 2£ 6 

x (2tt) + p 2 + p 3 - pa - p 5 - p 6 ) \ M 3 ^ 3 ? (6) 

where | ^3^3 | 2 is the squared amplitude for the 3-quark to 3-quark scattering. One 
consequence of such numerical simulations is to unfortunately break the locality of the 
transport equation and thus Lorentz covariance. Nevertheless, the particle subdivision 
technique [20,21] can restore the Lorentz covariance [11] by means of the transformation 

/->/' = if, 

I M 2 ^ 2 | 2 -| M'm | 2 =| M 2 ^ 2 I 2 /£, I | 2 ^| MU | 2 =| | 2 jt\ 

in combination with the replacement 

1/1/2(1 - / 3 )(i - h) - hh(i - h){\ - f 2 )] - hh - hh 

[/l/2/3(l - / 4 )(1 - / 5 )(1 - / 6 ) - hhhil - /l)(l - / 2 )(1 - / 3 )] - /1/2/S - /4/ 5 /6 

In principle, the Lorentz covariance is recovered while I goes to infinity. The practical 
value of i in codes equals 40. Six hundred and sixty-six quarks are generated from the 
initial distribution (5) and the average of 20 runs of Fortran codes with different random 
number sets is taken as a solution of Eq. (4). 

The solution of the transport equation at ^ so = 2 fm/c is shown by the dotted, 
dashed and dot-dashed curves in Fig. 4. The curves overlap and can thus be fitted to the 
Juttner distribution, 

/(P.*fao) = e ia/r_ A ( ? ) 

where the temperature of quark matter T = 0.59 GeV and fugacity A = 0.04. We get a 
thermalization time of t iso — t ini = 1.8 fm/c. For strongly interacting quark matter, the 
present value of the coupling constant o; s = 0.6 is allowed [25]. For weakly interacting 
quark matter, a s needs to be reduced. If the initial quark distribution and/or the coupling 
constant becomes smaller, the thermalization time of quark matter gets longer than the 
order of 1.8 fm/c. 

5. Summary 

We have studied the thermalization of quark matter including three-quark elastic 
scatterings in the transport equation. Squared amplitudes for triple-quark scatterings are 
derived at the tree level in QCD. Given an anisotropic quark distribution at the formation 
of quark matter, momentum isotropy can be shown by the transport equation solution. 
The history of quark matter's evolving into a thermal state can be described by the trans- 
port equation. The triple-quark scatterings give a variation of quark distribution function 



6 



comparable to the one that resulted from the quark-quark scatterings. Rapid thermal- 
ization can not be established by the triple-quark scatterings from such an anisotropic 
quark distribution which can appropriately describe quark matter initially created in cen- 
tral Au-Au collisions at the highest RHIC energy. Therefore, we also conclude that rapid 
thermalization of quark matter must depend on the gluon-gluon-quark, gluon-quark-quark 
and gluon-antiquark-quark scatterings which involve higher gluon number density than 
the quark number density and have larger squared amplitudes than the quark-quark- 
quark scatterings. Such very complicated works including the quark scatterings with 
gluons and/or antiquarks remain to be done. 
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Figure 1: Two-gluon-exchange induced scatterings of three quarks. 
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Figure 3: Quark distribution functions versus momentum in different directions while 
anisotropic quark matter is produced in the initial Au-Au collision. The dotted, dashed 
and dot-dashed curves correspond to the angles 6 = 0°,45°, 90°, respectively. 
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Figure 4: Quark distribution functions versus momentum in different directions while 
quark matter just arrives at thermal equilibrium. The dotted, dashed and dot-dashed 
curves correspond to the angles 6 = 0°,45°, 90°, respectively. The solid curve represents 
the thermal distribution function. 
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Table 1: For the entries we use the symbols (a u • • -,a 23 ) =(1024, 512, 2560, 1536, 2048, 
3072, 128, 256, 384, 768, 640, 896, 1152, 1280, 2304, 1792, 3328, 4096, 64, 448, 320, 192, 
832)/3456. 
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Table 1. Continued. 
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